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Photoproduction of D* and Jets at H1 



Zlatka Staykova on the behalf of H1 Collaboration 

Photoproduction events containing a charmed meson D** and two jets were investigated with 

the HI detector using the HERA II data sample. The D* meson was reconstructed in the decay 

channel, D +± — > — > K^n ± n ± . Jets were reconstructed using the inclusive k, algorithm and 

were selected if they have transverse momenta /? r (jet) > 3.5 GeV. One of the jets was associated 

with the D* meson itself, such that the jet originating from the parent charmed quark as the meson 

can be tagged. The phase space of the measurement is limited within central rapidity for the D* 

meson and the D* ., |tj| < 1.5 while the second jet was measured within, — 1.5 < TJ < 2.9. Single 
jet 

differential cross sections and double differential distributions were measured and compared to 
Leading Order Monte Carlo (MC) event generators, Pythia and CASCADE and with the Next- 
to-Leading order MC generator MC @NLO 
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D* and Jets at HI in yp 



1. Introduction and Motivation 

The dominant production process of charmed quarks in DIS is the boson gluon fusion process 
(BGF) (Fig. |I|), where a gluon splits into quark anti-quark pair and then interacts with a virtual 
photon. Such charm production is highly sensitive to the gluon content of the proton. Moreover 
charm mass m c ~ 1.5 GeV provides the hard scale for perturbative calculations in the regime of 
photoproduction where the photon virtuality approaches zero. In summary charm quarks in photo- 
production provide good testing ground for pQCD calculations. 

Higher order processes are often approximated by 
parton showers, where the successive emitted partons e — e /— 



from the proton side are ordered in a given kinematic 
quantity. Two major approaches are used, the DGLAP 
and the CCFM evolution equations. The basic differ- 
ence between them is the quantity used for the order- 
ing, in the case of DGLAP the transverse momentum 
of the partons is used, while in CCFM the emission 




angle with respect to the incoming gluon. 

Figure 1: A Feynman diagram of the bo- 
In photoproduction the photon is quasi-real can son gluon fusion process, 
split into partons. These events are called resolved 

photoproduction and are treated in terms of photon parton density functions (pdf) and evolution 
equations can be also applied to the photons. This is of particular need for the DGLAP based mod- 
els [JTJl- The process when the photon interacts as a pointlike object is called direct photoproduction. 

A previous measurement ^ of photoproduction of D* mesons and two jets at HI, was found 
to be highly sensitive to the two different parton shower approaches. In the present measurement a 
larger data sample is used and the phase space was extended towards larger rapidity for the second 
jet. 

Monte Carlo Models 

The presented measurement was compared to two leading order Monte Carlo (MC) programs, 
Pythia ^ and Cascade [Q]. In Pythia the parton showers are implemented according to the 
DGLAP evolution equations. It generates also resolved photon processes. Two different modes of 
Pythia were used: in the first case the matrix elements were calculated explicitly for heavy quarks, 
and the final prediction is a mixture of two different samples, direct and resolved processes. For 
that case the SAS 2D LO photon pdfs were used as well as the CTEQ 6M NLO proton pdfs. In 
the second mode the generator calculates the matrix elements for massless quarks. In this case the 
photon pdfs were GRV-G LO and the CTEQ 6L LO proton pdfs were also used. The CASCADE 
MC generated the parton showers according to the CCFM evolution equation with the set AO unin- 
tegrated proton pdfs. A next-to-leading order MC generator, (MC@NLO) [§] was also compared 
to the data. It is a full NLO matrix elements for heavy quark photoproduction matched with par- 
tons showers. The parton showers are implemented according to the DGLAP evolution equations, 
and were taken from the Herwig MC program and the CTEQ 6.6 proton pdfs were used. The 
uncertainty of the calculation was estimated varying the factorisation and renormalisation scales. 
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2. Experimental Setup, Event Selection and Event Reconstruction 
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Figure 2: The mass difference Am 

m{Kn) -m(Knn s i 0W ) 



In this measurement the HERA II data sample col- 
lected with the HI detector ^ was used, correspond- 
ing to an integrated luminosity of 5£ = 93.4pb _1 . The 
measurement was performed in untagged photopro- 
duction. The scattered electron escapes detection and 
the event kinematic variables were reconstructed using 
the hadronic final state (HFS). For the online trigger- 
ing of the events the fast track trigger (FTT) [^] was 
used. Charged particles are reconstructed online and 
offline and combined into a D* meson candidates using 
the mass hypothesis. The D* was reconstructed in the 
channel: D* ± D ^ K^i^ri^. The D* mesons 
were selected if they have transverse momentum p t > 2.1 GeV limited in the central rapidity range 
by the central tracking device |t]| < 1.5. The jets were reconstructed with the inclusive k t algorithm 
[]|] in the laboratory frame in the energy recombination scheme. The D* was treated as a leading 
particle, which means that the four vectors of the decay products of the meson were replaced by the 

four vector of the D* itself in the HFS definition. Such, the D* . could be identified. The second 

jet 

hardest jet in the event is denoted as the other jet. The jets were selected if they have transverse 
jet 

momentum of p\ > 3.5GeV in the central rapidity range for the D* , \r\{D* )| < 1.5 to be con- 
jet jet 

sistent with the D* selection. The other jet was measured in the range — 1.5 < r\ (Other jet) < 2.9. 
Finally a cut on the invariant mass of the jets My > 6 GeV was applied. 

The number of signal events was determined from a fit to the D* mass difference distribution 
Am = m(Kn) — m{Knn), (Fig. ||). For the signal the asymmetric Crystal Ball ^ function was 
used. The fit was performed in each bin of the measurement. Within the selection about 4000 D* 
mesons were found. 

Various sources of systematic uncertainty were investigated leading to total systematic uncer- 
tainty of about 10%. The dominant sources are the luminosity measurement, track finding and the 
trigger efficiency uncertainties. 



3. Results of the Measurement 

The cross sections are measured in the phase space as summarised above. The results are 
presented in Fig.|| - Fig.||. The cross sections are represented by the black markers, where the 
inner error bars shows the systematic uncertainty and the outer error bars are the total uncertainty. 

The kinematic quantities of the D* and the jets (not shown here) are well reproduced by the MC 
generators. In Fig. [3] the differential cross section of D* and two jets in photoproduction are shown 
as a function of the azimuthal angle difference between the two jets Aq> (top left); the invariant mass 
of the remnant in the event Mx (top right); the average transverse momentum of the di-jet pair pft 
(bottom left) and the longitudinal momentum fraction of the photon carried by the jets x r (bottom 
right). Here, the shape of the azimuthal angle difference A(p is not well reproduced by three models. 
The back-to-back region A(p close to 180°, where the jets are balanced in p t and contribution from 
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Figure 3: The differential cross sections for D* and two jet production in photoproduction as a function of 
Acp between the jets, Mx, py and Xy. 



higher order gluon radiation is not expected, is well modeled by Pythia massive. In the region 
towards small Aq>, further radiation is expected from momentum conservation, the data are well 
predicted by CASCADE. The same can be observed in the average transverse momentum of the 

jet pair pi . The invariant mass of the remnant M x = (p + J— (v\ + V2)) 2 , where v i / 2 are the four 
vectors of the two jets, is found to be well described by the models. The x r distribution is well 
described by the models with an exception of the lowest bin, where the contributions from events 
containing resolved photons is dominant. 
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Figure 4: The Acp distribution for two bins of Xy, compared to LO MC programs 



The Acp distribution was also measured for two bins of x y (Fig. Q). For the direct photon 
case, Xy > 0.75, within one sigma the distribution is fairly well modelled by CASCADE. It can be 
seen that at small Acp the model is slightly above the data. At the smallest Acp bin both PYTHIA 
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models are underestimating the data by a factor of 5. In the resolved case, x y < 0.75, CASCADE 
significantly underestimates the data in the back-to-back region, with a difference of factor of 4 
while at the tail of the distribution towards small A<p the model describes the data. The predictions 
of Pythia for the smallest Aq> bin is somewhat better than for the direct case. The data were 
also compared to MC@NLO. It was found that for the low x y < 0.75 (not shown here) region 
the predictions are below the data but for the direct photon x 7 > 0.75 case the shape is very well 
reproduced (see Fig. ||). 
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Figure 5: The A<p distribution for the high Xy region compared to MC@NLO 



4. Conclusions 

Cross sections of charm photoproduction with jets are highly sensitive to different parton 
shower models. A new measurement of D* and two jets in photoproduction with the HI detector 
was presented and differential cross sections and double differential distributions were measured 
with a larger data sample and in an extended phase space than a previous measurement by HI. 
Two new observables Mx and pj were introduced and measured. The measured cross sections 
were compared with predictions using different MC generators based on different parton shower 
models. It was found that the cross sections as a function of Aq> shows very different shape for the 
data and the models. 
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